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We propose two statistical models for description the metal-insulator phase transition coupled 
with paramagnetic-ferromagnetic phase transition in manganites of the type La\- x Sr x MnOz. The 
first one based on the competition of small polarons and delocalized carriers. In the second one the 
conductivity appears as a result of overlapping of large polarons. We compare our models with the 
experimental data. 
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The interplay of the metal-insulator (M-I) and 
ferromagnetic-paramagnetic (FM-PM) transitions in 
mixed metal oxides is one of the most important and 
least understood topics in material physics. The early 
explanation of the combined metal-insulator (M-I) and 
ferromagnetic-paramagnetic (FM-PM) transitions was 
based on the Double- Exchange (DE) mechanism [Q. In 
this mechanism the exchange interaction between spins is 
provided by conduction electrons. The electron tunneling 
amplitude from one site to another depends strongly on 
orientation of spins on these sites. Therefore, magnetic 
phase transition causes simultaneous drop in resistivity. 

Recently, Millis et al. , |ij have shown that the Dou- 
ble Exchange itself does not explain the transition from 
the conducting to insulating phase. They argued that 
Jahn- Teller distortions are strong enough and lead to for- 
mation of polarons. The high temperature phase in their 
picture contains presumably polaronic carriers and is in- 
sulating jjj. In the low temperature phase delocalized 
carriers appear together with the magnetization. This 
idea is supported by experimental observations of large 
Jahn- Teller distortions It was confirmed experimen- 
tally that small polarons are responsible for transport at 
high temperature [^). For reviews on the experimental 
state of art see M. 

Quantitative description of the phenomenon was based 
on a phenomenological Hamiltonian proposed by Millis, 
Schraiman and Mueller (MSM) H which describes elec- 
trons coupled via elastic deformations and interacting 
with localized spins 

H = ]T ^r+a^rMr + a)-£ J H Sr ■ ^(r)a^(r) 

r,a r 

+ E (f - 9Qr ■ ^(r)rV(r) - /^(r^r)) (1) 
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Here t Tir+a is the nearest neighbor tunneling ampli- 
tude, ?/>(r), i/^(r) are fermion operators, Q r are local 
Jahn- Teller-displacements, a and r are spin and pseudo- 
spin Pauli matrices, Jh is the Hund energy, g is the 



electron-phonon coupling constant and k is the elastic 
coefficient. Jh is assumed to be the largest energy scale 
in the problem. In this approximation electron spins fol- 
low adiabatically the direction of localized spins on a site 
The authors || treated their model Hamiltonian in 
the mean field approximation. For the number of carriers 
per site n = 1 they have found an insulator phase with 
strong Jahn- Teller distortion. For n < 1 the high temper- 
ature phase in the frameworks of this approximation re- 
mains conducting, though the conductivity can decrease 
rapidly above transition temperature. In particular, no 
remarkable magneto-resistance have been found at ex- 
perimentally optimal value n = 0.3. These results allow 
a following interpretation. The fermion variables can be 
eliminated explicitly from the partition function based 
on the Hamiltonian (0). After elimination we arrive at 
an effective free energy for variables Q and S: 

f = -TV In (l + exp (t - J - Q + £)) - - £ Q 2 (2) 

r 

where operator /t is diagonal 47V x 47V matrix in vari- 
ables r, a, r, whereas operator t is off-diagonal only with 
respect to variable r with non-zero matrix elements 
t J2 a <5r,r'+a; other two operators J and Q are off-diagonal 
with respect to variables a and r with matrix elements 
Jf/S r • a a ^i and gQ r ■ t t , t /, respectively. The mean-field 
(Landau) counterpart of the free energy (|^) reads: 

F L = aS 2 + b{S 2 ) 2 + AQ 2 + B(Q 2 ) 2 + cS 2 Q 2 . (3) 

This is well-known problem of two interacting order pa- 
rameters (see, e.g. 0). Its typical phase diagram con- 
tains 4 phases with S and Q either equal to zero or non- 
zero. It is tempting to identify phases with Q ^ as 
insulating and those with Q = as metallic. Then 4 men- 
tioned phases would be paramagnetic metal (PM), para- 
magnetic insulator (PI), ferromagnetic metal (FM) and 
ferromagnetic insulator (FI). However, this model be- 
ing applied for description of manganites with CMR has 
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two substantial failures: i) its high-temperature phase is 
metallic, whereas in real manganites it is insulating; ii) 
the transition from the PI to the FM phase proceeds at 
one point only, whereas in real systems it proceeds along 
a line. Therefore, the mean-field approximation applied 
by MSM from our point of view is insufficient to explain 
the real phase diagram in manganites. It happens be- 
cause any solution with the homogeneous displacement 
field Q does not describe formation of a localized objects 
- polarons. The fact that MSM obtained large magneto- 
resistance at n w 1 is readily explained as a result of 
splitting of the double-degenerate upper level by strong 
Jahn- Teller effect. But this effect does not work properly 
at smaller values of n, since the lower band is not fully 
occupied. 

In this article we propose simple statistical models in 
which polaron states are incorporated phenomenologi- 
cally. These models have the two most important phases 
found in the experiment: paramagnetic insulator (PI) 
and ferromagnetic metal (FM). In the first model which 
we call Small Polaron (SP) Model we start with the free 
energy written in the form: 

T= (-(l + m 2 ) + ^m 4 )tn c + Tm 2 +T LG (4) 

where Tlg = —n p s + T(n c lnn c + n p lnn p ) is the small 
polaron lattice gas free energy. n c is the density of con- 
ductance electrons, n p = n — n c is the density of lo- 
calized electrons (small polarons), n is the total elec- 
tron density; t is the bandwidth, — e is the small po- 
laron energy, m is the magnetization; the constant (3 al- 
lows to vary the saturation magnetization. In the insu- 
lating (polaron) phase m = 0. Minimizing free energy 
Eq.[I| with respect to m one finds equation for x — n c /n: 
xA + 9 ^2 x -i + Tin = 0. Depending on parameters, 
this equation has one, two or no roots. For A > and 
T = this equation has no roots. It means that the po- 
laronic phase is stable. For A < the FM phase is stable 
at T = 0. Since the polaronic phase always wins at large 
temperature, it implies that for A < the PI-FM phase 
transition takes place at any value of n.. The phase di- 
agram was calculated numerically. For (3 < l/(nln2) 
the transition temperature grows monotonously with n. 
If (3 > l/(nln2) the ferromagnetic phase exists not at 
any n, but only at the carrier density exceeding a criti- 
cal value n c < n. An interesting feature of this model is 
the possibility of re-entrant polaronic phase at low tem- 
perature. This peculiarity may be related to the reen- 
trant behavior of resistivity found in the experiment || . 
The energy parameters of the model e and t are in the 
range leV, but numerical values of the transition tem- 
perature are relatively small: ant < T c < bnt, where 
a » 0.01, b w 0.3. 

Though the SP model discussed above displays PI and 
FM phases and the colossal magneto-resistance, some its 



features seem unrealistic. It implies the first order phase 
transition with rather large discontinuities of the values 
m, n c and entropy for a reasonable set of parameters, 
whereas there is no experimental indication of remarkable 
hysteresis near the transition line. These discontinuities 
may be smeared out by a small number of impurities or 
other defects. We simulated the action of defects by aver- 
aging of the resulting values m(T/t,e),n c (T/t 7 e) over a 
small range (3 4%) of the parameter e. The results are 
shown in Fig.[l] for zero and small finite magnetic field. 
In this calculation we identified n^ 1 with the resistance. 
In reality they differ by a factor (mobility) which varies 
slowly, at least near the transition line. Another unrealis- 
tic feature of the SP model is too large scale of magnetic 
fields (this feature is shared with the MFA calculations 
of the MSM work). If one accepts the intersite tunneling 
constant equal to 0.6eV Jl|, the essential magnetic field 
at room temperature is about 10-100 times more than ex- 
perimental values (field in Fig.|l| corresponds to ~ 50T 
for t = O.GeV). 
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FIG. 1. Resistivity p(arbitrary units) vs T/t without and 
with (dashed line) effective magnetic field ho (see text). 

Our second model was inspired by the neutron scatter- 
ing experiments , [jl0| which showed the occurrence of 
large magnetic polarons (LP) coexisting with the small 
polarons (SP). This phenomenon lies beyond the scope 
of the simple SP model. Below we discuss an alternative 
mechanism of the PI-FM phase transition based on the 
LP. 

We employ Varma theory ||ll|| modified to incorporate 
the Jahn- Teller effect. The free energy of an individual 
carrier localized in a sphere of radius R reads: 

F(R) = ~W+-L-CTR 3 -^ (5) 

Here W oc Ms the half- width of the conductivity band, R 
is measured in the lattice constants, t and e ~ e have the 
same meaning as before; C is a numerical constant which 
can be roughly estimated as (4/3)7rln4 = 5.8. The latter 
figure is probably overestimated since it is found for the 
case of full saturation of magnetization in large polaron 



2 



which is scarcely realized. Further we assume C to be 
of the order of one. The first three terms in the r.-h. s. 
of eqn. ||were discussed in (yj. The last term describes 
the polaronic effect. It contributes a term proportional to 
|*(R)| 4 to the energy density where \&(R) is the electron 
wave- function (see, for example The contribution 

to the total energy is inverse proportional to the volume 



of the localized state. Exactly as in the work |11|, the 
localization on large scale of R originates from the en- 
tropy term TR 3 . The Jahn- Teller term can be neglected 
if t and Eq have the same order of magnitude and T is 
small. Minimizing F(R), we find the radius of the sta- 
ble or metastable polaron R p w (2t/3T) 1 / 5 . This result 
shows that the linear size of the LP changes very slowly 
with temperature. At reasonable values of t ~ Q.6eV jl3] 
and T ~ 0.01 -r 0.03ey it does not exceed 2.5. Thus, the 
large polaron is not so large. Still it contains between 30 
and 60 sites and its magnetic moment is large. 

The ratio of numbers of the LP and SP is given by 
the Boltzmann factor exp[(e — F(R p ))/T] where e is, as 
above, the energy of the small polaron. We introduce the 
value A p = e + W. If A p is positive, this ratio tends to 
a small number exp(— R p ) at high temperature. At some 
low enough temperature T p the LP density becomes equal 
to the inverse volume of the LP. This temperature can be 

found from equation: jfi- — 1.96 f J + In (nV p ) = 

where V p = (4/3)irR p is the volume of the LP. At this 
(or slightly higher) temperature large polarons overlap 
forming an infinite cluster. The percolation leads to 
a rapid decrease of resistivity. Indeed, the magnetiza- 
tion inside the cluster is homogeneous. The condition 
at which the spontaneous magnetization appears was 
shown to be tn c > T. The concentration of delocal- 
ized carriers n c in the cluster is equal to V^ -1 . Thus, 
tn c /T p « (i/T p ) 2 / 5 > 1. 

Normally the percolation transition proceeds contin- 
uously. The Coulomb repulsion prevents formation of 
macroscopic magnetic droplets before percolation. The 
Coulomb energy e 2 /eR must be compared with the ki- 
netic energy t/R 2 . Even for static dielectric constant 
e = 25 the Coulomb energy prevails. We expect that the 
screening is even weaker at small distances. 

It is important to note that the existence of magne- 
tization inside the cluster does not guarantee the total 
non-zero magnetization. The percolation and the sponta- 
neous magnetization appearance proceed simultaneously 
in 2 dimensions. The reason is that two infinite clusters 
for different spin orientations would unavoidably inter- 
sect. In 3 dimensions they can interpenetrates without 
intersection as long as their branches are thin enough. 
The spontaneous magnetization appears when the cross 
section of each infinite cluster increases. This fact was 
established for Ising model Q], but its geometrical ori- 
gin is quite general. In the frameworks of the LP model 
the Mott transition proceeds, generally speaking, at a 



temperature higher than the Curie point. 

An approximate solution of equation for T p is: T p ~ 
t(A p /2tf/ 2 . At A p « 0.5* we find T p ~ 0.03* » 200.fi:. 
This estimate looks reasonable. The logarithmic term in 
equation for T p provides a weak dependence of T p on the 
dopant concentration n. 

Statistics of LP can be constructed via a kind of the 
lattice model with the sites either occupied by a LP or 
empty (the elementary cell of this auxiliary lattice has the 
volume V p ). Introducing occupation numbers r r = 0,1 
and polaron spins S r 
Hamiltonian reads: 



5n r on sites, the proposed lattice 



*r+a ' r ' r+a 



Hi p = Jn r • n r 

r.a 

+ ^ AV r T r+a ~ ^2( h ' "r^r + fJ>T r ), (6) 



where J is the effective exchange constant, h = gfigHS 
is an effective magnetic field for polarons, K stays for the 
Coulomb interaction between neighboring polarons and 
fx = W — CW z ^T 2 / b is the polaron free energy. The 
mean-field counterpart of the Hamiltonian (j^j reads: 



Fi 



ip 



Jm 2 x 2 — TiS^r(m) — hmx — J-i(x) (7) 



where m — |(n r }| and x = (r r ), J = zJ, K = zK, 
z is coordination number. J-i{x) = Kx 2 — T[x\ux + 
(1 — x) In (1 — x)] — fix. Su{rn) is the spin entropy per 
site for the Heisenberg system. For large S we are inter- 
ested in, Snim) is determined by a following equations: 
S H {m) = In 25 - ln(l - (l-m)((m)) - (m + l)C(m), 
where ^(m) is determined by equation: coth£— 1/( = m. 
We calculated magnetization by minimization of free en- 
ergy (0) and resistance employing the effective medium 
approximation [IB] ]. The results for magnetization and 
large polaron density are shown in Fig.^. 



FIG. 2. Magnetization (line) and large polaron concen- 
tration (+) vs T/W for the set of parameters J = .2W 
K = 0.6W. 
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Despite of visible similarity of the LP mean- field model 
to corresponding SP model, the former displays the sec- 
ond order transition in a broad range of parameters. The 
resistivity show magnetic field dependence qualitatively 
similar to that in Figjlj However, due to large value of S 
for a polaron, the scale of magnetic field displaying the 
CMR is about 100 times lower than for the SP model and, 
thus, it fits satisfactory the experimental data. Transi- 
tion temperature in Fig.|^ is in qualitative agreement with 
estimate from the simple large polaron model of Eq||. 
The phase diagram of the model Eq.f^ displays region of 
the second order phase transitions for 0.5VL < K < 1.5W 
and for J < 0.2-^0.3IT with T p oc J. The Mean Field Ap- 
proximation (MFA) is insensitive to percolation thresh- 
old near which the phase transition occur. However, one 
can estimate correction to the MFA phase diagram by 
using results known from dilute magnetic models [fl6|| . 
Corresponding corrections decrease the transition tem- 
perature in the region with LP concentration larger than 
percolation threshold by « 10% 20%. More detailed 
description of the phase diagram will be given elsewhere. 

In conclusion, we proposed statistical models explic- 
itly incorporating polaron formation in manganites. One 
of them (the SP model) displays first order phase tran- 
sition from paramagnetic state with localized electrons 
to the ferromagnetic state with metallic conductivity in 
perovskites, the colossal magneto-resistance and the re- 
entrant resistivity in a range of variables. The second 
LP model is based on percolation of large polarons. It 
predicts a continuous isolator-conductor transition at a 
temperature higher than the magnetic transition. Ex- 
perimental data on neutron scattering, the magnetic 
and resistive measurements near the transition line al- 
low to give some preference to the LP model, though 
the further experimental and theoretical work is neces- 
sary. When this work was completed and presented at 
the Argonne Symposium on Poorly Conducting Oxides 
(Argonne, June 1998), R.H. Heffner kindly informed us 
that he and his colleagues were driven to a similar con- 
clusion on the existence of dynamic conducting clusters 
treating results of their experiments on muon spin relax- 
ation in La\- x C a x MnOz |Tij| . 
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